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Abstract Genome size changes significantly among 
taxonomic levels, and this variation is often related to 
the patterns shaped by the phylogeny, life histories and 
ecological factors. However, there are mixed evidences on 
the main factors affecting molecular evolution in animals. 
In this study, we used phylogenetic comparative analysis 
to investigate the evolutionary rate of genome size and the 
relationships between genome size and life histories (ie. 
hatchling mass, clutch size, clutches per year, age at sexual 
maturity, lifespan and body mass) among 199 squamata 
species. Our results showed that the evolutionary rate of 
genome size in Lacertilia was significantly faster than 
Serpentes. Moreover, we also found that larger species 
showed larger hatchling mass, more clutches per year 
and clutch size and longer lifespan. However, genome size 
was negatively associated with clutch size and clutches 
per year, but not associated with body mass we looked at. 
The findings suggest that larger species do not possess the 
evolution of large genomes in squamata. 
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1. Introduction 


Genome sizes vary considerably across taxa in organisms 
(Cavalier-Smith, 1978; Bennett and Leitch, 2005; Lynch and 
Walsh, 2007). This can be driven by the stochastic genetic and/ 
or genomic processes associated with spontaneous deletions 
and/or insertions, polyploidization, prolonged tandem repeats 
length, transposable elements number and genetic drift, but can 
be also shaped by natural selection (Ogata et al, 1996; Petrov, 
2001; Sun et al., 2012; Lynch, 2011; Whitney and Garland, 2010). 
In particular, genome size variations are mainly explained 
by two important mechanisms including the duplication 
events and the proliferations of noncoding elements (Neiman 
et al., 2015). Establishing the association between genome size 
variation and organismal complexity has puzzled many 
evolutionary biologists and as such remains a classic problem 
in biology (Gregory, 20052). 

Previous studies across taxa have revealed positive 
associations between genome size and cell size, nucleus 
size, developmental time, nutrient requirements, tissue 
differentiation, life cycle complexity and body size (Vinogradov, 
1997; Olmo and Morescalchi, 1978; Gregory, 2005a; Gregory, 
2001; Gregory and Johnston, 2008; Guignard et al, 2016). These 
positive associations have been suggested to be consequences of 
both the cytoplasm from more efficient mRNA transport and 
larger cells necessitating larger genomes based on structural 
causes (Cavalier-Smith, 1985). Smaller cells for instance, usually 
divide faster and have a higher metabolic rate, evidenced by 
a negative correlation between metabolic rates and DNA 
amounts in mammals and birds (Hughes and Hughes, 1995; 
Vinogradov, 1995; Gregory, 2002a; Hughes and Piontkivska, 
2005). However, a potential correlation which still needs to be 
explored between cell volume and genome size is body size 
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(Gregory et al. 2000). 

Body size variation is often determined either by cell size 
and cell number, or both combination in organisms (Hessen et 
al, 2013; Kozłowski et al, 2003). For plants and animals, genome 
size displays a positive association with cell size (Gregory et al., 
2000; Bennett, 1987; Gregory, 2005b). In addition, this correlation 
can be often linked to ecological factors. For example, genome 
size exhibits a positive association with body size in some 
invertebrates (e.g, amphipods, copepods, crustaceans) due to low 
metabolic rate and temperatures in cold waters (Rees et al., 2008; 
Angilletta et al., 2004; Timofeev, 2001; Jeffery et al, 2016; Leinaas 
et al, 2016). Genome size variation in frogs is indirectly affected 
by temperature and humidity as a result of its influence on the 
time of premetamorphic development (Liedtke et al., 2018). In 
birds, mitochondrial and nuclear of substitution rate in coding 
sequences reveal weak negative associations between the ratio 
of nonsynonymous and synonymous substitution rate and 
age at sexual maturity, lifespan and body mass associated with 
environmental factors (Weber et al., 2014; Lanfear et al., 2010; 
Nabholz et al, 2013), but it is not always the case (Figuet et al., 
2017). 

Squamata constitutes the class of vertebrates with a small 
genome size due to a lower fraction of transposable elements 
and shorter introns (Organ et al., 2007). However, this may be a 
misconception caused by overlooking GC-rich regions, which 
are often hard to access (Botero-Castro, 2017). Ploidy variations 
does not provide a major power of variation in genome size, 
and the whole-genome duplication events is not reported 
during the amniote evolutionary process in squamata (Van de 
Peer et al., 2009). Squamata mainly consists of two suborders 
(eg, Lacertilia and Serpentes) and displays complex life histories 
with prolonged developmental periods (hatching time), which 
likely constrains the variation of genome size because of a 
negative correlation between genome size and development 
time in invertebrates (Wyngaard et al., 2005). 

To examine the selective mechanisms underlying genome 
size variation in squamata, we first estimated the evolutionary 
rates of genome size between Lacertilia and Serpentes in 
squamata. We also expanded our extent to which genome 
size can be considered as a determinant of life histories by 
investigating the relationships between variation in genome 
size and life histories among 199 squamata species. We tested 
whether larger bodies can promote evolution of larger genomes. 


2. Materials and Methods 


2.1. Data collection The genome size of 199 squamata 
species was collected from genome size database (http://www. 
genomesize.com) (supplementary Information: Table S1). 
We extracted data on genome size for squamata species for 
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which information on life histories can be found (see below), 
and obtain their average C-value. We used average values of 
genome size when more than one measurement per species 
was available. To avoid possible errors due to several methods 
being used to quantify genome size (Hardie et al., 2002), we used 
parallel analyses on a subset of genome size. We confirmed 
species names using the NCBI taxonomy database, and 
collapsed/pruned all synonyms from the phylogenetic tree. We 
rebuilt the phylogenetic tree using time-calibrated molecular 
phylogeny by Pyron et al. (2013) (Figure 1) and examined 
difference in the evolutionary rate of genome size between 
Lacertilia and Serpentes. Finally, we compiled information on 
hatching time, hatchling mass, clutch size per year, clutch size 
and body mass (see details in De Smet, 1981; Feldman et al., 2016; 
Allen et al, 2017) and age at sexual maturity, lifespan from the 
AnAge databases (https://genomics.senescence.info/species/) 
(supplementary Information: Table SI). 


22. Statistical analyses The complementary approaches were 
used to evaluate the evolution rate of genome size for three 
suborders. For each suborder, we assessed phylogenetic signal 
using the phylosig function in the package of pli ytools in RStudio 
v.312 (Revell, 2012). We then used the Blomberg's K (Blomberg 
et al., 2003) in which genome size variation comparing on a 
null model is assumed genome size evolution under Brownian 
motion (BM) model. We also used the Pagel's A (Pagel, 1999) 
in which phylogenetic signal is estimated on the basis of the 
phylogenetic dependence of genome size. K - 1 indicated 
genome size evolved as expected under a BM model, while K 
»lorK < 1 indicated less or more phylogenetic signal than 
expected under a BM model, respectively. We used Blomberg’s 
K and Pagel's À to estimate the phylogenetic signal and found 
qualitatively similar results (Table S2). 

We used the fitContinuous function in the R package-Geiger 
(Harmon et al, 2008) to compare genome size evolution on the 
basis of Brownian motion, Ornstein-Uhlenbeck and Early-burst 
models between the two suborders. Following the suggestions 
by Simmons and Fitzpatrick (2016), BM model of genome 
size evolution was regarded to be the best model due to small 
sample size. Moreover, to compare differences in evolutionary 
rate of genome size between the two suborders, we modified 
a likelihood method where a phylogeny can directly compare 
on the Brownian evolutionary rate (7^) of genome size (Adams, 
2013). 

To examine associations between genome size and life 
histories, we used the phylogenetic generalized least squares 
models where the phylogenetic structure of the model 
residuals was considered in the caper package (Orme et al, 2012; 
Huang et al., 2020). We used phylogenetic scaling parameter 
À to estimate the phylogenetic influence on the associations 
between genome size and life histories based on a maximum- 
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likelihood approach (Pagel, 1999). The scale of A-values ranges 
from zero (i.e, phylogenetic independence) to one (i.e, complete 
phylogenetic non-independence) (Freckleton et al., 2002; Mai 
et al, 2019). We log-transformed life histories to linearize 
associations and used the phylogenetic tree of squamata species 
to correct for phylogenetic dependence (Mai et al., 2020). To 
test the associations between body mass and life histories, we 
treated body mass as response variable, hatchling mass, clutches 
per year, clutch size, age at sexual maturity and lifespan 
as predictor variables using the multivariate phylogenetic 
generalized least squares. To test whether genome size exhibited 
a association with body mass, we treated body mass as predictor 
variable, genome size as response variable, and hatchling mass, 
clutches per year, clutch size and lifespan as covariates using 


the multivariate phylogenetic generalized least squares. 
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3. Results 


The average value of genome size was 2.11 pg, ranging from 
119 to 393 pg among 199 species of squamata. Genome size in 
Lacertilia tended to be larger than that in Serpentes (Figure 2). 
The evolutionary rate of genome size in Lacertilia was faster 
than that in Serpentes (Table S3). 

The multivariate phylogenetic generalized least squares 
model indicated that body mass was positively associated with 
hatchling mass, clutches per year, clutch size and lifespan 
among 199 species of squamata (Table I) . The genome size was 
not associated with body mass when the effects of hatchling 
mass, clutches per year, clutch size and lifespan were removed 
(Table 2). We also found negative correlations between genome 
size and clutch size or clutches per year (Table 2). 

For Serpentes in particular, body mass was positively and 
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Figure 1 The phylogenetic tree of the 199 species of squamata used in the comparative analysis. 
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significantly associated with hatchling mass and clutch size, but 
not with clutches per year, age at sexual maturity and lifespan 
using the multivariate phylogenetic generalized least squares 
model (Table S4). However, there was no association between 
genome size and body mass when removing the hatchling 
mass and clutch size effects (Table S5). For Lacertilia, body mass 
was significantly associated with hatchling mass, clutches per 
year, clutch size and lifespan (Table S4). When the influences 
of hatchling mass, clutches per year, clutch size and lifespan 
were removed, we found no association between genome size 
variation and body mass (Table S5). 


Table 1 The associations between body mass and life histories across 
199 species of squamata. Phylogenetic scaling parameters (superscripts 
following A denote P-values of likelihood ratio tests against models with 4 
- 0 and À = 1, respectively). 


; Body mass 
Predictors 1 a : r2 P 
Hatchling mass «0.001^?'5 1.001 9.543 0.791 «0.001 
Clutch size 0.811 5.059 0.516 «0.001 
Clutches per year 0.530 2.293 0.180 0.031 


Age at sexual maturity -0.372 -1.541 0.090 0.136 
Lifespan 0.481 2.947 0.266 0.007 


Table 2 The associations between genome size and life histories across 
199 species of squamata. Phylogenetic scaling parameters (superscripts 
following 4 denote P-values of likelihood ratio tests against models with 2 
- 0 and 4 =1, respectively). 


Genome size 


Predictors 1 A ; E P 

Body mass «0.0017?! 0.090 1.865 0.110 0.073 
Hatchling mass -0.087 -1.392 0.065 0.175 
Clutch size -0.231 -3.920 0.354 «0.001 
Clutches per year -0.146 -2.124 0.139 0.043 
Lifespan 0.028 0.565 0.011 0.576 


4. Discussion 


Our results showed that genome size evolution in Lacertilia 
evolved significantly faster than that in Serpentes among 199 
species of squamata. We found positive correlations between 
body mass and hatchling mass, clutches per year, clutch size, 
and lifespan. However, genome size was not associated with 
body mass when correcting for the effects of part life histories. 
For Lacertilia and Serpentes, genome size did not show a 
association with body mass. 

Differences in transposable element accumulation rates 
in animals experienced may lead to substantial variation 
in genome size among species (Chalopin et al, 2015; Gibbs 
et al., 2004). For example, a number of DNA obtained by 
transposable element accumulation with strong changes among 
lineages, are counteracted by loss of DNA on the basis of large 
segmental deletion in birds (Kapusta et al, 2017). For 199 species 
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Figure 2 Genome size difference between Lacertilian and Serpentes for 
199 species of squamata. 


of squamata, the rate of transposable element accumulation 
can also explain the marked variation in genome size, ranging 
from 1.19 to 3.93 pg. 

The evolutionary history of genome size in amphibians 
has been one of gradual, time-dependent variation (Brownian 
motion; Liedtke et al., 2018). In this study, evolutionary model- 
fitting showed that genomes in Lacertilia and Serpentes evolved 
under a shared processes of Brownian motion. The common 
ancestor of extant squamata was predicted to have similar size 
in genome in Lacertilia and Serpentes. We inferred that genome 
size in squamata evolved gradually as a function of time 
(Brownian motion). Herein we found that the evolutionary 
rate of genome size in Lacertilia evolved faster than Serpentes. 
Palaeontological data and genomic evidence display a similar 
pattern (Pyron et al, 2013). 

There are evidences that phylogeny is likely to promote the 
influences of genome duplications and transposons on genome 
size evolution in animals (mammals: Tang et al., 2019; insects: 
Alfsnes et al., 2017). For example, genome size is phylogeny- 
dependent when 4 » 09 in all life-history traits is reported in 
mammals (Tang et al., 2019). However, phylogeny displays a 
weak correlation with genome size in crustaceans (Alfsnes et al., 
2017). Likewise, there is a weak association between genome size 
and phylogeny among 240 species of birds when A < 0.564 is 
recorded in all life-history traits (Yu et al., 2020). We found that 
genome size was not associated with phylogeny, suggesting that 
the phylogeny did not a strong power in driving transposons 
and duplications of genome in squamata. 

Genome size variation can be explained by the more 
mechanistic and/or short period effects which is regarded as 
the proximate causes. Moreover, the evolutionary powers (ie, 
selection), regarding as the ultimate causes, can also explain the 
genome size variation (Hessen et al, 2013; Alfsnes et al., 2017; Yu 
et al., 2020). For birds, variations in genome size are positively 
related to the length of developmental period (Kapusta et al., 
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2017; Yu et al., 2020), providing evidence for the associations 
between life histories and genome size evolution. Indeed, 
genome size displays markedly and directly effects on cell size 
and cell replication rate (Gregory, 2002b), so larger genomes 
are expected to be positively correlated with larger egg size and 
smaller clutch size. However, large datasets have indicated that 
variations in genome size are not associated with offspring 
number and size in mammals (Tang et al., 2019) and life history 
complexity of amphibians (Liedtke et al., 2018). In this study, 
there were negative correlations between genome size variation 
and life histories such as clutch size and clutches per year in 
squamata, suggesting that less offspring number or larger 
offspring size can promote evolution of larger genomes. 

Body mass is positively associated with genome size in 
vertebrates (Liedtke et al., 2018; Tang et al, 2019; Yu et al, 2020) 
and invertebrates (Gregory et al., 2000; McLaren et al., 1989; 
Hessen and Persson, 2009; Alfsnes et al., 2017). Such positive 
associations between cell size and genome size (Gregory, 2005a; 
McLaren and Marcogliese, 1983) have indicated that variations 
in body size among the related species can partly respond to 
variation in cell size (Hessen et al., 2013). Indeed, genome size 
exhibits positively correlations with body mass in birds and 
mammals (Tang et al., 2019; Yu et al, 2020). Across 199 species 
of squamata, there were not associations between genome size 
and body mass, suggesting that diversity in genome size was 
not response of variation in cell size. 

In conclusion, we illustrated the relationships between 
genome size and life histories in squamata. The hatching time, 
hatchling mass, clutch size per year and clutch size cannot 
shaped the genome size variation, and species with larger 
bodies did not possess larger genomes in squamata. Our future 
research would need more species to reveal the relationships 
between genome size evolution and life histories. 
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Appendix 


Table S1 TSpecies, body mass (g), genome size (pg), hatchling mass (g), clutch size, clutches per year among 199 species of squamata from the references of 
De Smet (1981), Feldman et al. (2016), Allen et al. (2017), and age at sexual maturity (years) and lifespan (years) from AnAge (https://genomics.senescence.info/ 
species/). 


Species Suborder Genome size Body mass Hatchling mass Clutch size Clutches per year Ageatsexualmaturity Lifespan 
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(Continued Table S1) 


Species Suborder Genome size Body mass Hatchling mass Clutch size Clutches per year Ageatsexualmaturity Lifespan 
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(Continued Table S1) 


Species Suborder ^ Genomesize Body mass Hatchling mass Clutch size Clutches per year Ageatsexualmaturity Lifespan 


Mabuya mabouya Lacertilia 127 12.82 0.69 410 
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(Continued Table S1) 


Species Suborder Genome size Body mass Hatchling mass Clutch size Clutches per year Ageatsexualmaturity Lifespan 
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Xenochrophis piscator Serpentes 168 13160 082 3450 100 NA 900 
Xerotyphlops vermicularis Serpentes 196 545 NA 550 NA NA 620 
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Table S2 Evaluation of phylogenetic signal in genome size examined. 


Suborder Blomberg's K Pagel’s A 

K P-value À InLA InLA = 0 (P-value) InLA = 1 (P-value) 
Lacertilia 0.418 0.001 0.606 131.30 103.480 (<0.001) 111.327 (<0.001) 
Serpentes 0.228 0.235 <0.001 69.814 69.814 (1.000) 59.969 (<0.001) 


Note: All analyses were conducted using log,)-transformed data. Phylogenetic signal was estimated by using Blomberg’s K and Pagel's 4. K and 
À of 1 are indicative of genome size evolution following a Brownian motion model. For Blomberg's K, significant P-values show phylogenetic 
signal in the genome size. For Pagel’s 2, the maximum likelihood estimate of à was compared using likelihood ratio test to models were å = 0 
and = 1, respectively. Significant P-values show if À differs significantly from models assuming no phylogenetic signal (A = 0) or strong phy- 
logenetic signal (4 = 1). 


Table S3 Comparison of model parameters and fit for each suborder examined under Brownian motion, Ornstein-Uhlenbeck and Early-burst 
evolutionary models. 


Suborde Brownian motion model Ornstein-Uhlenbeck model Early burst model 

iid c InL AICc A; 0 a Int  AlCc A 00 a InL AICc A 0 
Lacertilia 20* 1133 -21856 2363 740? 001 12419 -24219 000 100 -100% 11133 -21646 2573 259? 
Serpentes 30* 5997 -11576 2021 409? 006 7116 -13597 000 100 -100^ -5997  -11358 2239 137? 


Note: The Brownian rate parameter, 7", and rate of evolutionary change parameter, a, are presented for the BM, and EB models, respectively. 
For all models the maximum likelihood estimates (InL) and sample size corrected Akaike Information Criterion (AICc) values are presented. 
To compare model fits, for each suborder we report the value of delta AICc, A, and the Akaike weights, w, which indicate the strength of 
evidence for each model. 


Table S4 Associations between body mass and life histories for the two suborders in squamata using phylogenetic generalized least squares models. 
Phylogenetic scaling parameters (superscripts following 2 denote P-values of likelihood ratio tests against models with A = 0 and A = 1, respectively). 


Body mass 
Suborder Predictors z 2 
À [2 t R P 
Hatchling mass «0.001 799?! 0.666 4.630 0.573 «0.001 
Clutch size 0.979 2.999 0.360 0.009 
Serpentes Clutches per year 0.229 0.640 0.025 0.531 
Age at sexual maturity 0.355 0.623 0.024 0.542 
Lifespan 0.416 1.222 0.085 0.239 
Hatchling mass «0.001'7^?! 1.186 15.999 0.850 «0.001 
Clutch size 0.978 7.001 0.522 «0.001 
Lacertilia Clutches per year 0.854 4.923 0.139 «0.001 
Age at sexual maturity -0.016 -0.079 0.085 0.937 
Lifespan 0.415 2.953 0.017 0.005 


Table S5 Associations between genome size and life histories in squamata using phylogenetic generalized least squares models. Phylogenetic scaling 
parameters (superscripts following A denote P-values of likelihood ratio tests against models with A = 0 and A = 1, respectively). 


Suborder Predictors i Z Cenom m R P 
Body mass «0001? 0201 0.469 0.006 0.642 

Serpentes Hatchling mass —0.038 -0.791 0.017 0434 
Clutch size -0.012 -0.155 <0.001 0.878 
Body mass 02460 0.021 0510 0.005 0.612 
Hatchling mass -0.020 -0.359 0,002 0721 

Lacertilia Clutch size -0.054 -0.864 0.013 0.391 
Clutches per year —0.045 —0.680 0.009 0.500 


Lifespan —0.006 -0.153 <0.001 0.879 


